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ABSTRACT 

Strongly-magnetized, rapidly-rotating neutron stars are contenders for the central engines of 
both long-duration gamma-ray bursts (LGRBs) and hydrogen-poor super-luminous super¬ 
novae (SLSNe-I). Models for typical (~ minute long) LGRBs invoke magnetars with high 
dipole magnetic fields (If \ > 10 15 G) and short spin-down times, while models for SLSNe-I 
invoke neutron stars with weaker fields and longer spin-down times of weeks. Here we iden¬ 
tify a transition region in the space of !f\ and birth period for which a magnetar can power 
both a long GRB and a luminous SN. In particular, we show that a 2 ms period magnetar 
with a spin-down time of ~ 10 4 s can explain the observations of both the ultra-long GRB 
111209 and its associated luminous SN2011kl. For magnetars with longer spin down times, 
we predict even longer duration (~ 10 6 s) GRBs and brighter supernovae, a correlation that 
extends to Swift J2058+05 (commonly interpreted as a tidal disruption event). We further 
show that previous estimates of the maximum rotational energy of a proto-magnetar were too 
conservative and energies up to E nrdx ~ 1 - 2 x 10 53 ergs are possible. The magnetar model 
can therefore comfortably accommodate the extreme energy requirements recently posed by 
the most luminous supernova ASASSN-151h. The high ionization flux from a pulsar wind 
nebula powering ASASSN-151h may lead to an ionization break-out X-ray burst over the 
coming months, which would be accompanied by an abrupt change in the optical spectrum. 
We conclude by briefly contrasting millisecond magnetar and black hole models for SLSNe 
and ultra-long GRBs. 
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1 INTRODUCTION 

In addition to canonical long-duration gamma-ray bursts (LGRBs) 
with characterstic durations of tens of seconds, a handful of bursts 
lasting longer than 1000 seconds have been discovered jLevan et al.| 
|2014| l. These “ultra-long” GRBs (ULGRBs) share many properties 
with LGRBs of normal duration, including rapid variability, large 
isotropic radiated energies close to a solar rest mass, and an asso¬ 
ciation with star-forming galaxies at cosmological distances (see 
|Levan||2015[ for a review). It remains debated whether ULGRBs 
are simply the longest lasting members of a single, continuous 
LGRB population ( [Virgili et al.|2013[[Zhang et al.|20 1 A\ or whether 
they represent a distinct class with potentially different progenitors 
( |Boer et al.|2015[|Levan et al.|2014| >. 

Most LGRBs are accompanied by stripped envelope, hyper- 
energetic supernovae (SNe) jWoosley & Bloom|2006j ), making it 
natural to ask whether the same is true for ULGRBs. The redshift 
of GRB 121027A was too high to search for a SN. while the af¬ 
terglow of GRB 130925A was heavily extinguished by dust. The 
afterglows of GRB 101225A and 111209A both showed blue to 
red color evolution several weeks after the burst, consistent with a 
SN brightening above the power-law afterglow {Levan et al.|2014[ 


Greiner et al. 2015| l. However, the spectrum of these thermal tran¬ 
sients extended further into the ultraviolet than other GRB SNe. 
The lack of clear spectral features was explained as Doppler line 
blending resulting from high ejecta velocities of > 20,000 km s~' 
( |Greiner et al.|20l5l >. 

If ULGRBs result from the core collapse of massive stars then, 
as with other LGRBs, central engine models fall into two classes: 
a hyper-accreting stellar-mass black hole (BH) l [Woosley||1993| > 
or a neutron star (NS) with a millisecond rotation period and an 
ultra-strong magnetic field (a “millisecond proto-magnetar”; |Usov| 
11992[ >. In magnetar models the GRB jet is powered by the rotational 
energy extracted by electromagnetic torques, and hence the GRB is 
generally weaker and longer in duration for magnetars with weaker 
magnetic fields. A BH engine, by contrast, is powered only so long 
as it continues to accrete stellar debris. The durations of ULGRBs 
challenge BH models because the majority of the bound debris re¬ 
turns on the gravitational free-fall time of the stellar envelope (e.g., 
Perna et al.| |2014| l. This is usually only minutes for the compact 
Wolf-Rayet progenitors of normal LGRBs, but can be months or 
longer for blue super-giants (|Quataert & Kasen|2012||Woosley &| 
|Heger||20l2| |Gendre et al.|20 1 3| . However, the lack of hydrogen 
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spectral features in the SNe which accompany ULGRBs may dis¬ 
favor such extended stellar progenitors [Greiner et al .|2015) . 

A long-lived central engine is also a leading explanation for 
the hydrogen-poor class of “super-luminous supernovae” (SLSNe 
Type I, or SLSNe-I), stellar explosions with peak luminosities 10 
to 100 times higher than normal core collapse SNe ( jQuimby et ah] 
2011; Gal-Yam 20121 and which cannot be powered by the radioac¬ 
tive decay of 56 Ni. In engine-powered models for SLSNe-I, the cen¬ 
tral compact object is again either a millisecond proto-magnetar 
(Kasen & Bildsten 2010||Woosley|2010] > or an accreting BH jDex-| 
|ter & K asen 2013 1 . The long-lived nature of the central engine, 
rather than just its total energy budget, is critical to enhancing the 
SN luminosity: energy released too early is degraded by adiabatic 
expansion of the ejecta before being radiated. A connection be¬ 
tween SLSNe-I and LGRBs is suggested by similarities in their 
host galaxy properties ( |Lunnan et al.|2 014; Leloudas et al.|2015] >. 

The discovery that SN 201 lkl following GRB 111209A was 
significantly over-luminous as compared to normal GRB-SNe pro¬ 
vides a potential link between ULGRBs and SLSNe-I, leading 
Greiner et al. < [2015| ) to suggest a magnetar origin for this event. 
This interpretation, however, also raises many questions: If a mag¬ 
netar (or an accreting BH) can separately produce both an LGRB 
and a SLSNe-I, can one produce both simultaneously? Why are 
the SNe that accompany ULGRBs different than the optically- 
selected SLSNe-I? Why are GRB-SNe usually not as luminous as 
SN 201 lkl? 


Adding fuel to the debate was the recent discovery of 
ASASSN-151h, a SLSNe-I with the highest peak luminosity (> 


2 x 10 45 ergs s _1 ) of any SN known to date. Its total radiated energy 
of ss 10 52 erg represents a significant fraction of the maximum rota¬ 
tional energy of a NS, posing a challenge to the magnetar model for 
SLSNe-I (and, indeed, to any model). In a similar way, the gamma- 
ray and afterglow kinetic energies of the brightest LGRBs detected 


by Fermi are inferred to approach or exceed a few 10 5 - ergs (Cenko 


|et al.|20n) . This apparent fine tuning of the spin rate so close to 
its maximum has generally been used to argue against the magne¬ 
tar model for LGRBs, although the clustering of GRB-SNe kinetic 
energies near the same scale has been used to argue in favor of a 
magnetar engine (Mazzali et al.|2014). 


In light of GRB 111209A/SN201 lkl and ASASSN-151h, as 
well as the rapidly growing sample of SLSNe-I (e.g.,|Tnserra et al.| 
]20T3[ [Nlcholl et al.||2014> , we revisit the range of observational 
signatures of magnetar birth, placing GRBs and SLSNe-I within a 
common framework. We show that “transition” events like LGRB 
111209A/SN2011kl are a natural consequence of the model, pro¬ 
vided that a relativistic jet can escape the star. Furthermore, al¬ 
though the energy of ASASSN-151h was extreme, we show that 
previous estimates of the maximum energy of a magnetar have been 
underestimated. 


2 ENERGETICS OF PROTO-MAGNETAR SPINDOWN 


The rotational energy of a NS of gravitational mass M ns is given by 

£ rot - IQ 2 /2 x 2.5 x 10 52 (M n JL4M o ) 3/2 P m 2 erg, (1) 

where P = 2n/Q. = P ms ms is the rotational period and / = 1.3 x 
10 45 (M ns /1.4M Q ) 3/2 gem 2 is the NS moment of inertia (Lattimer 
& Schutz 2005 ). The maximum allowed rotational energy, £ max , 
corresponds to the minimum spin period, P min , set by the mass- 
shedding limit. In ^4] we show that P m j n typically varies from x 1 
ms for M x 1.4M 0 to 0.7 ms for M ~ 2M 0 (depending on the NS 
equation of state), such that £ max > 10 53 erg for massive NSs. 


The pulsar loses rotational energy at a rate which, for an 
aligned force-free wind, is given by (Contopoulos et al.| 1999[) 


T s d — 


fi 2 Q 4 


—- 1.7 x 10 s 


R- P - ' 


l + — 

Ld 


ergs 


( 2 ) 


where fj. = B d RjJ, s * s th e dipole moment, fid = 10 15 Bi 5 G is the 
surface equatorial dipole fieltj]] fi NS = 12 km is the NS radius, and 



1=0 


147s(M/1.4M 0 ) 3/2 B7 5 2 fi 2 ms 


(3) 


is the initial spin-down time. The spin-down luminosity exceeds the 
value given by equation 0 for several seconds after core bounce 
due to the neutrino-heated wind ^Thompson et al.]2004j [Metzger] 
jet al.|20lTj>, but this correction is minor for ULGRBs. 


2.1 Powering Gamma-Ray Bursts 

The magnetar outflow, though approximately isotropic on small ra¬ 
dial scales, is collimated into a narrow jet along the polar axis via 
the confinement provided by the surrounding stellar envelope jL'z-| 
jdensky & MacFadyen||2006| Bucciantini et al.||2007| l. This colli- 
mation may be facilitated by the anisotropic stress induced by the 
strong toroidal magnetic field in the nebula separating the mag¬ 
netar wind from the surrounding star (Bucciantini et al.|2009) . It 
remains uncertain whether such a strong, ordered magnetic field is 
maintained in the face of non-axisymmetric instabilities |Begelman| 
[T998l|Porth et al,|2013||Mosta et al.|20T4> . 

The formation of a stable jet does not guarantee that the jet 
will escape the star over the timescale f sd of peak spin-down 
power. A jet of luminosity Lj and half-opening angle 0j requires a 
time fesc ~ M e jV e jC0 2 /Lj to escape the stellar ejecta of mass M e j and 
velocity v e j. The condition that f esc < f sd for Lj « L sd | f=0 ~ E m Jt si 
translates into a maximum opening angle for jet escape, 

*®) 1/2 ’ jetescape (4) 

( jQuataert & Kasen|2012j > where the second equality assumes that 
E mt x M e jV 2 /2, as applies if the kinetic energy of the explosion is 
dominated by the magnetar rotational energy. Condition 0 is not 
easily verified by current simulations, hut it still provides a consis¬ 
tency check on the opening angle derived by modeling ( §3.1[ >. 

The power of a jet that cleanly escapes the star will approxi¬ 
mately equal the peak spin-down power of the NS, i.e. Lj =s L sd |, =0 
(eq. j2j) [Bucciantini et aL| |2007j >. The observed isotropic y-ray 
luminosity is then Li so « e y ff l L s a\ l= o, where f b x fjr/ 2 is the beam¬ 
ing fraction and e y is the radiative efficiency. The LGRB duration is 
most naturally associated with the spin-down time over which L sd 
is roughly constant, i.e. T y x f sd (eq. (3)), assuming the jet propa¬ 
gation time through the star is negligible in comparison. 


2.2 Powering Supernova Light Curves 


The SN explosion ejects a mass of M e j with an initial kinetic energy 
Eox = 10 51 erg. As the ejecta expands, its thermal energy E evolves 
as a function of time t: 


dE 

dt 


“~Vej - Lsn + Lnj + L sd , 
^ej 


(5) 


1 Note that our definition of B d, which agrees with that most commonly 
used in t he puls ar community, is lower than that value adopted by |Kasen~&] 


Bildsten 


(2010 by a factor of Vl2 = 3.46. 
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Figure 1. Parameter space of dipole magnetic field B& and birth spin period Pq, calculated for M e \ = 3M©, A/ni = 0.2M 0 , k = 0.2 cm 2 g -1 , E^o = 10 51 erg 
and M = 1.5M©. Dashed blue contours show the peak spin-down luminosity (top panel) and the characteristic spin-down timescale (bottom panel), the later 
potentially associated with the LGRB duration. Solid black contours show the peak SN luminosity Lsn (top panel) and peak timescale ?sn (bottom panel). 
Superluminous (SLSNe), powered by NS rotational power, have Lsn > 10 43 erg s -1 . Dark green circles show the NS parameters inferred by previous fits to 
other SLSNe-I by|Nicholl et al. (2014) and |Chatzopoulos et al.|(2013) (we caution that the best fit ejecta mass was not usually equal to the 3 M Q assumed in 
calculating the contour lines). Red circles show fits from this paper to GRB 111209A/SN201 lkl (Greiner et al.|2015) , ASASSN-151h (Dong et al.|2015) , and 
Sw J2058+05 (Pasham et al.|2015) . GRB111209A/SN201 lkl is constrained by the requirements: SN peak luminosity Lsn = 3 x 10 43 ergs -1 (solid black line), 
SN peak time fsN = 15 days (dashed black line), ejecta velocity v e j « 2 x 10 4 km s -1 (solid green line), and a GRB of duration T y = 4d ~ 1.4 x 10 4 s (dashed 
blue line). ASASSN-151h is constrained by: SN peak luminosity Lsn = 3 x 10 45 erg s -1 and SN peak time ?sn ~ 35 days. 
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where the first term accounts for PdV losses, R e j and v e j = dR^/dt 
are the mean radius and velocity of the ejecta, respectively. The 
kinetic energy E k = M c ,jV 2 ./2, increases as the result of PdV work, 
dE k /dt = (£7_R c j)v e j, as required by energy conservation. 

The term L SN = E/t A is the radiated luminosity, where f d = 
(3/rMej)/(4TCv e ji) is the photon diffusion timescale, k is the optical 
opacity, and L Ni ~ 1.3 x lO^ergs” 1 (M Ni /0.2M o )exp(-f/8.8d) is 
the heating rate due to 56 Ni decay, where M Ni is the 56 Ni mass. The 
latter varies from ~ 0.05 - 0.1M o in normal Type Ib/c SNe 
to M Ni « 0.3 - 0.5 M g in GRB SNe ( [Nomoto et al,|[2003^ . The 
opacity of k = 0.1 - 0.2 cm 2 g _1 is set by electron scattering and 
Doppler-broadened lines. The last term in equation § accounts for 
energy input from the NS (eq. (2j), which we assume thermalizes 
with unity efficiency (at late times t » f sd ; see ©• 

3 THE DIVERSITY OF TRANSIENTS FROM 
MAGNETAR BIRTH 

We survey the landscape of GRB/SN transients produced by mil¬ 
lisecond NS birth for a range of dipole fields S d ~ 10 12 -10 16 G and 
initial spin periods Po — 0.7 - 10 ms. All solutions are calculated 
for M c j = 3M 0 , M Ni = 0.2M Q , k = 0.2 cm 2 g _1 , and M = 1.5M 0 . 

Figure |T] overviews the 6 d - P 0 plane, with the top and bot¬ 
tom panels showing contours of GRB/SN luminosities and peak 
timescales, respectively. Dashed blue lines show the maximum jet 
luminosity Lj = L sd |, =0 (top panel; eq. |2j ) and the initial spin-down 
time f sd (bottom panel; eq. Magnetars in the far upper-left cor¬ 
ner of the diagram produce high jet luminosities and short spin- 
down times t si < 100 s, characteristic of the (beaming-corrected) 
luminosities and durations of normal LGRBs. Moving to longer 
birth periods and lower B d , the spin-down time increases and the 
jet becomes weaker. If the jet is stable and can escape the star for 
t si > 10 3 - 10 4 s, such events could produce ULGRBs. 

Solid black lines show the peak SN luminosity, Lsn (top 
panel), while dashed lines show the time of the SN peak, t SN (bot¬ 
tom panel). For both the highest magnetic fields of B d > 10 15 G, 
as well as very low B d /large Po, which likely characterize the birth 
of ordinary pulsars, the SN light curve is predominantly powered 
by M 'Ni with a characteristic duration of f SN s» 30 days. However, 
for P < 10 ms and intermediate values of i? d ~ 10 12 - 10 l3 G the 
SN luminosity is enhanced to a superluminous level by energy in¬ 
put from the pulsar ( |Kasen & Bildsten|201Ofr . The “sweet spot” at 
B d js 10 13 G and Po ~ 1 ms corresponds to a NS containing as 
much rotational energy as possible, with a spin-down time f sd com¬ 
parable to the ejecta diffusion time at peak light. Green points show 
the best-fit magnetar parameters for other SLSNe-I (JNicholl et al.| 
|2014||Chatzopoulos et al.|2013), while red points show the fits from 
this paper for GRB111209A/SN2011kl, ASASSN-151h, and Swift 
J2058+05. 

Figure [7] shows that the region of normal duration LGRBs 
overlaps with the Ni-powered SN regime, consistent with the stan¬ 
dard interpretation that GRB-SNe are powered by radioactive Ni 
(as in other Type Ic SNe; e.g., Nomoto et al. 2003). However, 
ULGRBs occupy a region of parameter space that overlaps that 
of SLSNe-I, especially for magnetars with short birth periods of 
Pa < 2 ms. Some of the longest GRBs could thus be accompanied 
by very luminous SNe l |Quataert & Kasen|2012| l. 

3.1 GRB 111209A/SN 2011kl 

SN 201 lkl peaked at a luminosity of L p k » 3 x 10 43 erg s _1 on a 
timescale of 14 days ([Greiner et al .[2015| ). This fast rise, coupled 
with the high ejecta velocity (> 20,000 km s _1 ) required to explain 



10 1 10 2 10 3 10 4 10 5 10 6 10 7 


tsd = T r (s) 

Figure 2. Peak SN luminosity, Lsn. as a function of the magnetar spin-down 
time t sd (eq. 0) for a range of models (black circles) with B d e [10 12 ,10 16 ] 
G and Pq e [1,10] ms. For short spin-down times, rotational energy is in¬ 
jected too early to significantly increase the SN luminosity above its mini¬ 
mum value set by 56 Ni heating (which can vary moderately from the value 
of Mni = 0.2M q assumed here). However, for longer f sd the maximum 
allowed value of Lsn increases. Assuming that the prompt GRB duration 
approximately equals the spin-down timescale, i.e. T y f sd , we show for 
comparison normal LGRB-SNe (blue circles; |Cano|2013) and the ULGRBs 
(red circles) 11 1209AjGreiner^af J2015) , 101225A [Levan et al.|2014I , 
and Swift J2058+05 jPasham et al.|201 57 


the spectrum, requires a low ejecta mass of a few solar masses. Fig- 
ure|T]shows that GRB 111209A and SN 201 lkl are simultaneously 
explained by a magnetar with B d j p 


3 x 10 14 G and P * 2 ms 


for AL c j x; 3M g . The preferred model has two free parameters (for 
fixed Mej) and lies at the intersection of four constraints: the peak 
duration and luminosity of the SN (dashed and solid black lines), 
the mean velocity of the ejecta (green line), and the GRB duration 
of T y » f sd « 1.4 x 10 4 s (dashed blue line). 

The model then predicts a beaming-corrected jet luminosity of 
Lj » 10 48 erg s _1 (solid blue line). The observed av erage isotropic 
gamma-ray luminosity of L y ~ 4 x 10 50 


2015) then implies a jet opening angle of 6 S « (2Ljg ) ,/L r ) 1,/ » 
0.05, assuming a typical value of e y = 0.5 ( jZhang et al.||2007| >. 
This is consistent with the maximum opening angle of (9j, max ~ 0.17 
for which the jet is capable of escaping the star (eq. (4), assuming 
v e j ~ 0.06 c). If the ejecta mass had been closer to the average 
value of « 10 M g inferred for SLSNe-I ( |Nicholl et al.|2015] >, then 
the jet opening angle for GRB 111209A would be much closer to 
the critical value. That a successful jet emerged from this explosion 
might therefore be the result of its lower than average ejecta mass. 

Figures [T] shows that ULGRBs from NSs with weaker mag¬ 
netic fields or longer spin periods than that responsible for GRB 
111209A should also produce dimmer SNe. As shown in Figure[2] 
the maximum peak SN luminosity should correlate with the burst 
duration for large values of the latter. Indeed, the SN associated 
with the ULGRB 101225A of shorter duration ( T y » 7 x 10 3 s) had 
a peak absolute magnitude M B « -19 ( [Levan et al.|2014| , closer to 
normal GRB SNe. Also following this trend is the extremely long 
X-ray transient Swift J2058+05 ( T y ~ 10 6 s; Pasham et ah[2 015 i, 
which, although generally interpreted as a tidal disruption event, 
could also be a core collapse event ( [Quataert & Kasen||20 1 2[ . A 
weaker correlation is expected for low values of T y due to varia¬ 
tions in the baseline luminosity set by the amount of 56 Ni produced 
in the explosion. However, it cannot be discounted that some nor- 
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mal LGRB-SNe could be at least partially engine-powered, low¬ 
ering the required 56 Ni mass. The engine’s contribution to the SN 
would also be underestimated by our model if the spin-down power 
is higher at early times than predicted by eq. 0 due, for instance, 
to a time-dependent NS magnetic dipole (or, in the case of BH- 
powered events, a complex accretion history). 

3.2 ASASSN-151h 


ASASSN-151h had a peak luminosity of Lsn > 2 x 10 42 erg s * 1 

10 52 ergs (Dong et al. 2015). 


and a total radiated energy of 
Figure [2] shows our best-fit “magnetar” model for ASASS-141h of 
B a = 1.2 x 10 13 and P 0 = 1.2 ms, calculated again assuming an 
ejecta mass of M c j = 3 M e . The solution is constrained by fitting 
the peak luminosity and timescale of the SN, resulting in an ejecta 
velocity of = 2.4 x 10 8 cm s _1 . It is challenging to power such a 
luminous event without also placing a comparable amount of ki¬ 
netic energy into the ejecta via PdV work. A large kinetic energy 
for ASASSN-151h is suggested by the presence of a broad feature 
in the optical spectrum at 4200 A with a width of 250 A, corre¬ 
sponding to a velocity spread of v e j ~ 2 x 10 4 km s _1 . 


4 MAXIMUM NEUTRON STAR ROTATIONAL ENERGY 

Once accounting for both radiation and kinetic energy imparted to 
the ejecta, ASASSN-151h requires an energy of > 2x 10 52 erg, close 
to the maximum rotational energy of a 1.4 M a NS. Such fine tuning 
appears to push the NS model to its limit by requiring ASASSN- 
151h to occupy a ‘pin prick’ in B d - P 0 parameter space. ASASSN- 
151h-like events are, however, very rare: |Dong et al.| ( |2015) ) esti¬ 
mates a rate 10-100 times lower than average SLSNe-I, which are 
already only ~ 10 -4 of all core collapse SNe (Quimby et al. 2011). 
Furthermore, previous estimates of the maximum rotational energy 
are too conservative when more massive NSs are considered. 

Figure [3] shows the maximum allowed NS rotational energy 
Emax (black line) as a function of the gravitational mass, calculated 
for a series of solid-body rotating NSs using the rns code (Ster- 
Igioulas & Friedman|1995| l. We assume a parametrized piecewise 
polytropic EOS, as described in |Margalit et al.||2015| >, with an adi¬ 
abatic index of T = 3 above the break density of p\ = 10 14 ' 7 g cm 3 
and pressure P i = 3.2 x 10 34 dyn cnr 2 . The particular EOS used 
in Fig. [3] results in a radius of 10.6 km and gravitational mass of 
3Tmax(0 = 0) ~ 2.24M 0 . Solid body rotation is established within 
minutes or less after NS formation due to the radial redistribution of 
angular momentum by magnetic stresses and gravitational waves. 

Figure [3] shows that a maximum rotational energy of E max » 
10 53 erg is reached for a NS with M x M mdx (Q. = 0). A solid- 
body rotating NS can possess a mass above this limit if it ro¬ 
tating sufficiently rapidly to be stabilized by centrifugal forces 
(a “supramassive” NS). Although the maximum rotational energy 
set by the mass-shedding limit continues to increase with M > 
M md x(f2 = 0) along the supramassive branch, the minimum rota¬ 
tional energy needed to stabilize the NS from collapse, 2t co u (red 
line), increases even faster with M. Hence, the maximum avail¬ 
able rotational energy, A£ max = E mdX - £ co n to be extracted de¬ 
creases for M > M mdX (Q = 0). The peak of AE mdX thus occurs at 
M » M mdX (Cl = 0). Although Figure [3] highlights our result for a 
particular EOS, we find that this maximum energy lies within the 
relatively narrow range of 0.9 - 1.65 x 10 53 erg across a wide range 
of EOS consistent with the lower limit set on the maximum mass 
of a non-rotating NS i jAntoniadis et aH2013| ). This maximum ro¬ 
tational energy could in principle be limited by gravitational wave 
emission for the highest rotation rates. However, we find that the 



M„ s (M 0 ) 


Figure 3. Maximum extractable rotational energy from a NS, A/T max = 
E mdx - Feu (black solid line), as a function of the NS mass M, where E mdX 
(blue dashed line) is the maximum rotational energy at the mass-shedding 
limit and F co u (red dotted line) is the minimum rotational energy required 
for support of a supramassive NS against collapse to a BH. The structure 
of the solid-body rotating NS is calculated using the rns code assuming a 
parametrized piecewise polytropic EOS with an adiabatic index T = 3 above 
the break density of p\ = 10 14 7 g cirT 3 at a pressure of P\ = 3.2 X10 34 dyn 
cnf 2 . We find that the maximum value of A E mdx lies within the relatively 
narrow range of 0.9-1.65x 10 53 erg across a wide range of I ’ P\ consistent 
with constraints on the maximum mass of a non-rotating NS. 


ratio of kinetic energy to gravitational binding energy, T/W, of the 
mass-shedding sequence generally remains below the critical value 
of (T/VT) max = 0.13 required for the growth of secular instabili¬ 
ties, driven by gravitational radiation reaction or viscosity ( |Lai &| 
Shapiro 1995). Thus, E mdx » 10 53 erg represents a new, relatively 
conservative upper limit on the energy permitted by the magnetar 
model, alleviating the tension produced by ASASSN-151h and the 
most energetic LGRBs. 


5 DISCUSSION AND CONCLUSIONS 


Although ASASSN-151h and GRB111209A/SN2011kl can be ex¬ 
plained by the magnetar model, we have made several implicit as¬ 
sumptions. A critical one is that the jet cleanly escapes the star 
at early times < t si , placing a large fraction of its energy into the 
GRB jet, only to become trapped at later times » t sA , such that 
the remaining rotational power thermalizes behind the ejecta with 
high enough efficiency to power the SN. Such a transition is possi¬ 
ble if the jet becomes less stable, or less effective at maintaining an 
open cavity through the expanding ejecta, as the spin-down power 
decreases from its peak initial value. 

It is useful to consider what conclusions can be drawn from 
the SLSNe-I population as a whole. Although many of the cur¬ 
rent sample of SLSNe-I fie far from the ULGRB region (Fig. [TJ, 
a few are sufficiently close to produce an ULGRB of duration 
T y > 10 4 s should the jet escape the star. Such high energy tran¬ 
sients could easily have been missed due to relativistic beaming 
away from our line of site, especially giv en th e beaming fraction in¬ 
ferred for GRB 111209A of ff ~ 400 (5 3.1 1 . Are the rates consis¬ 


tent? SLSNe-I occur at rates comparable to those of normal LGRBs 
( [Quimby et al.|20TT| . However, marginally over-luminous events 

like SN2011kl are more common than more luminous SLSNe- 

I, while ULGRB s may be intrinsically rare compared to normal 
GRBs, even after accounting for the bias against detecting lower 
luminosity events. The SLSNe-I SCP06F6 did show an X-ray out¬ 
burst of luminosity L x ~ 10 45 er g s 1 several months after the ex- 
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plosion ( |Levan et al.||2013) , but this may have signaled ionizing 
radiation escaping through the ejecta instead of a true jet (see be¬ 
low). Orphan afterglow emission from an off-axis jet would peak 
at radio wavelengths on a timescale of months after the explosion, 
motivating a systematic radio search following nearby SLSNe. 

Figure [I] also shows that the current SLSNe-I sample does 
not uniformly sample the NS parameter space capable of produc¬ 
ing SLSNe, with a clear bias present for events with higher val¬ 
ues of Bd. In principle, millisecond pulsars with fields as weak 
as B d ~ 10 12 G should produce SLSNe-I, albeit ones with slowly 
evolving light curves. This dearth of low-B d events may be intrinsic 
if NSs acquire strong magnetic fields as the result of a dynamo in 
the proto-NS phase caused by rapid rotation {Thompson & Duncan| 
|1993| >. It is also possible that assumptions of the SN model break 
down at low-B d , in particular regarding the efficiency with which 
the pulsar wind thermalizes its energy inside the SN ejecta. As 
the ejecta expands and its density decreases, non-thermal radiation 
from the pulsar wind nebula can ionize its way through the ejecta 
i jMetzger et al.|20L4| ), allowing pulsar power to escape directly be¬ 
fore thermalizing into optical radiation. High energy gamma-rays 
also more readily escape at late times due to the decreasing y - y 
and y — e cross sections at photon energies » m e c 2 . If the ther- 
malization efficiency is reduced for the low-B d pulsars, which pro¬ 
duce SLSNe-I that nominally peak at low Thomson optical depths 
T es < 30, this could help explain the apparent dearth of low-B d NSs 
in the SLSNe-I population. 

ASASSN-151h is the only SLSNe-I that we find lies below 
t cs( ? sn) < 30. Such luminous events can produce an “ionization 
break-out” signal as soon as a few months after the SN ([Metzger] 
|et al.|2014j >, characterized by an abrupt onset of soft X-ray emis¬ 
sion of luminosity L x ~ 10 42 - 10 44 erg s _1 and a drastic change 
in the optical spectrum to one characterized by nebular emission 
lines. Such a break-out event is not guaranteed to occur (especially 
if the ejecta mass is high), but we strongly encourage high-cadence 
X-ray and optical spectroscopic coverage of ASASSN-151h over 
the coming months, as its discovery would be a smoking gun for 
engine-powered models of SLSNe. Other than their unusually long 
durations, ULGRBs resemble the bulk of the LGRB population in 
many of their other properties. Therefore, if a single ULGRB can 
be unambiguously associated with magnetar birth, this strengthens 
the hypothesis that many (or all) LGRBs are magnetar-powered 
(Thompson & Duncan|1993| |Thompson et al.|2004||Metzger et al.| 
|201 l||Mazzali et al.|2014) 

Finally, much of the physics driving high energy transients and 
SLSNe from magnetar birth applies equally well if the engine is in¬ 
stead an accreting BH ( [Dexter & Kasen|20 1 3f , provided that the 
magnetar spin-down time is replaced by the period of peak accre¬ 
tion and the r 2 decay is replaced by the commonly expected r 5/3 
decay in the mass fail-back rate from the stellar envelope. Correla¬ 
tions between direct engine activity and SLSN luminosity such as 
those shown in Figure [2] should hold more generally. Perhaps the 
best discriminant between NS and BH models for ULGRBs is the 
prediction of hydrogen in the SN ejecta in the BH case due to the 
requirement of an extended progenitor star. Although hydrogen is 
not detected in the current sample of SNe accompanying ULGRBs, 
Doppler broadening of the lines could reduce their strength. Future 
observations and more detailed modeling are thus required to place 
more meaningful constraints. 

If the jet power tracks the magnetic flux accumulated by the 
BH instead of just its accretion rate ( [Tchekhovskoy & Giannios| 
|2015[ >, a large progenitor star with a long free-fall time might not 
necessarily be required to explain ULGRBs. In this case, LGRBs of 


shorter(longer) duration could result from stars with largerfsmaller) 
magnetic flux. However, in this case a seemingly fine-tuned corre¬ 
lation between the magnetic flux and the normalization of the fall 
back rate would be required to explain the similar total jet energies 
in normal and ULGRBs. 
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